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Abstract 

The aim of Dawn mission is the acquisition of data from orbits 
around two bodies, (4) Vesta and (1) Ceres, the two most massive 
asteroids. 

Due to the low thrust propulsion, Dawn will slowly cross and transit 
through ground-track resonances, where the perturbations on Dawn 
orbit may be significant. In this context, to safety go the Dawn mission 
from the approach orbit to the lowest science orbit, it is essential to 
know the properties of the crossed resonances. 

This paper analytically investigates the properties of the major 
ground-track resonances (1:1, 1:2, 2:3 and 3:2) appearing for Vesta 
orbiters: location of the equilibria, aperture of the resonances and pe- 
riod at the stable equilibria. We develop a general method using an 
averaged Hamiltonian formulation with a spherical harmonic approxi- 
mation of the gravity field. If the values of the gravity field coefficient 
change, our method stays correct and applicable. We also discuss the 
effect of one error on the C20 and C22 coefficients on the properties of 
the 1:1 resonance. These results are checked by numerical tests. We 
determine that the increase of the eccentricity appearing in the 2:3 
resonance is due to the J22 coefficient. Finally we numerically study 
the probability of the capture in resonance 1:1. Our paper reproduces, 
explains and completes the results of Tricarico and Sykes (2010). 

Keyword Vesta - Gravitational perturbations - Resonance - Space- 
craft operations - Analytical study 



1 Introduction 



After the explorations of the planet of our Solar System, the new space mis- 
sions are aimed at asteroids and comets, moving from fly-bys to rendezvous 
and orbiting. For these missions the space probe could have an operating 
orbital period close to the asteroid one. A challenge for the navigators of 
these missions is to predict the orbital environment around the asteroid and 
to derive pre- mission plans for the control of these orbits. 

Se veral authors have already been interested in these problems. Scheeresl 
(| 19941 ) studied the stability of a spacecr aft in synchronous orbits with the 



asteroid rotation. IScheeres et al.1 (119981) presented the dynamics of orbits 



close to the asteroid 4179 Toutatis. 



Rossi et al.l (|l999T ) numerically stud- 
ied the orbital evolution around irregular bodies with two methods. The 
first is based on a representation of the central body in terms of "mascons" , 
and the second model considers the central body as a polyh edron with a 



variable number of triangular faces. iHu and Scheeresl (|2002l ) investigated 
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the spacecraft motion about a slowly rotating asteroid (as compared to the 
space craft orbit perio d) wit h the second deg r ee an d order of the gravity 
field. ILara et alJ ifcoid ) and IRussell and Laral (l2009h were interested in the 



motio n of an orbiter around Enceladus. More recently, Tricarico and Svkesl 



<j20ld > studied the dynamical environment of Dawn around Vesta. 



The asteroid (4) Vesta is one of the biggest minor planets of the main 
belt. Originally discovered i n 1807 by Heinr i ch Pipers, it is one of the t argets 
of the space mission Dawn ( Ravman et al. . 20061 : Russell et al. . 20071 ). that 



will perform three different altitude science orbits before flying to Ceres. 
Dawn was successfully launched on September 27, 200 7 and will arrive to 
Vesta in July 2011. In the preparation of this mission, Tricarico and Svkesl 



(l2O10l ) aroposed a numerical survey of the orbits around Vesta and empha- 



sized particularly the problem of the 1:1 ground-track resonance. 

We revisit their study in adding a more complete numerical exploration 
and analytical explanations of the observed phenomena i.e. location, aper- 
ture and period of the stable equilibria of the ground-track resonances. The 
interest of an analytical investigation is the easiness and the fastness to 
refind the results when the data change. 

Our paper is organized as follows. In section [2] we briefly recall the 
specification of the Dawn mission and we determine the main forces acting 
on Dawn orbiting Vesta. Then we select the most important perturbations 
th at we use in our analytica l models. In section [3] we rebuild the results 
of Tricarico and Svkesl ( 2010l ) and we numerically localize the ground-track 



resonances (the commensurabilities between the orbital period of the probe 
and the rotational period of the asteroid); we measure the aperture of these 
resonances and we determine the period at the stable equilibrium of the 1:1 
resonance. In section [J] we recover the results of section [3] but in a purely 
analytical way. We develop a general method using an averaged Hamiltonian 
formulation with a spherical harmonic approximation of the gravity field. 
We also discuss the effect of an error on the C20 and C22 coefficients on the 
properties of the 1:1 resonance. In section [5] we show that the increase of the 
eccentricity in the 2:3 resonance is due to the J22 coefficient. In section [6] 
we study numerically the thrust and the probability of capture in the 1:1 
resonance and we detail the results of Tricarico about the trapping. Finally 
we summarize and draw our conclusions. 
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2 Main perturbations acting on Dawn around Vesta 



2.1 Dawn mission 

Dawn will be the first mission to orbit two main belt asteroids: Vesta and 
Ceres. They are the two most massive asteroids in the main belt and they 
are particularly interesting because they span the region of space between 
the rocky inner body Solar System and the wetter bodies of the outer So- 
lar System. They have survived intact from the earliest days of the Solar 
System. Hence they can provide the opportunity for acting as windows to 
the conditions present in the first few millions years of the Solar System 
existence. 



For a co mplete overvi e w of t he Dawn mission, we refer to lRayman et al. 



(l2006l l and iRussell et all The science observations will b 



e concen- 



trated in three campaigns around Vesta, each one conducted in a different 
circular, Sun-synchronous quasi-polar orbit. The first science orbit, Vesta 
Science Orbit 1 (VSOl), has a radial distance of 2 700 km (corresponding to 
an equatorial altitude of about 2 400 km) with a period of 58 hours. Upon 
completion of VSOl, the thrusting will resume for the 3- week transfer to 
Vesta Science Orbit 2 (VS02). The VS02 distance is 950 km (an equatorial 
altitude about 670 km), where the orbital period is 12 hours. This part of the 
mission is named HAMO, namely High Altitude Mapping Orbit. The trans- 
fer from VS02 to VS03 will take 30 days. At a radial distance of 460 km 
(equatorial altitude of about 180 km), VS03 is the lowest orbit planned and 
its orbital period is 4 hours. This last scientific orbit is named LAMO for 
Low Altitude Mapping Orbit. The transition between the different phases 
are made by a Solar electric propulsion with low thrust. 

2.2 Main perturbations 

The major perturbations acting on a probe orbiting an asteroid of the main 
belt are due to the discrepancies of the asteroid gravitational field, to the 
gravitational perturbations induced by the Sun, Mars and Jupiter as well as 
to the dir ect radiation pressure of the Sun. W e use the following simplified 
relations ( Montenbruck and Eberhard , 2000l ) to characterize the norms of 
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the accelerations induced by these 3 forces: 



„ j^n j 

O-nm = { n + 1) ^ ~^ V + ^nm ) 

= ^3b r (I) 

U 3b 

a rp = C r — Pq 
m 

where a nm is proportional to the coefficient of degree n and order m of the 
spherical harmonic development; a Sb corresponds to the acceleration induced 
by the third body (Sun, Mars or Jupiter) and a rp is the acceleration due 
to the direct radiation pressure. R e is the equatorial radius of the asteroid, 
Cnm and S nm are the normalized coefficients of the gravitational potential; // 
and fi represent the gravitational constants of the asteroid (Vesta or Ceres) 
and of the third body respectively. Pq = 4.56 x 10 N/m is the radiation 
pressure for an object located at a distance of 1 AU from the Sun and 
C r is the reflectivity coefficient appearing in the direct radiation pressure 
acceleration. A/m is the area-to-mass ratio, d 3b is the distance between 
the object and the third body and r is the distance between the probe 
and the central body (asteroid). Th e A/m of Dawn is approximately equal 
to 0.04 m 2 /kg ( Tricarico and Svkes . 2O10l ). Th e values of the coefficients 



Cnmi S nm that we use can be found in Table 1 of lTricarico and Svkesl (l2O10h 
as well as the values of = 2.36 UA and R e = 300 km. 

We draw the order of magnitude of the main perturbations (JTJ) as func- 
tions of the "planetocentric" satellite distance r in Fig. [Tj As a consequence 
of the location (~ 2.5 UA) and the shape of the asteroid (Vesta and Ceres) 
the order of magnitude of the perturbations is relatively different from the 
magnitude of the perturbations for a satellite of the Earth. For a space 
probe orbiting Ceres or Vesta, the perturbation due to the discrepancies of 
the gravitational field is several orders of magnitude greater than these due 
to the third body (or to the radiation pressure) perturbations. This remains 
true for large values of the probe semi-major axis. On the contrary, in the 
case of a stationary satellite of the Earth (GEO), the Luni-Solar perturba- 
tions are comparable in size to the Earth oblateness ones. 

Therefore, later in this paper, we shall only keep the perturbations due 
to the gravitational field of the asteroid. 
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Figure 1: The orders of magnitude of the main perturbations acting on a 
probe around Vesta and Ceres (left and central panel) and, comparatively, 
around the Earth. The dashed vertical lines show the "geo" stationary and 
"asterostationary" semi-major axis. 



3 Vesta dynamical environment: numerical study 



The effects due to a non-spherical asteroid shape are quite complex. Sev- 
eral a uthors gave differe nt solutions to modelize the asteroid gravitational 
field. iRossi et al l dlOQflh us ed the mascons or the polyhedrons approxima- 
tions. Scheeresh 19941 ) and Compere et al. ( 2011 ) used a tri- axial ellipsoid 
representation. Let us note that the spherical har monic representation of the 



;ep 
k3, 



1994; Pick et al 



gravi ty field is uniformly convergent for a < c\/2 ( Balmin 

where a and c are the biggest and smallest semi-axis of the asteroid. 
This is the case of Vesta. Then the major gravitational perturbations of the 
asteroid shape onto the orbit of an orbiter can usually be characterized by 
the spherical harmonics. For example, to characterize the main effect of the 

asteroid shape, it is convenient to use the primary oblaten ess term C20. 

Fu rthermore, using a spherical harmonic representation. iTricarico and Sykes 
(boioh ;ested several models for the mass density distribution of Vesta. They 
concluded that Dawn dynamics depends primarly on Vesta shape and only 
secondarily on the details of the interior structure. Therefore we decide to 
use the uniform mass density scenario of Tricarico with a spherical harmonic 
development of the gravity field. 
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As previously stated, the dynamics of a spacecraft orbiting within ~ 
1 000 km from Vesta is dominated by the gravitational field of Vesta. As a 
consequence the system of differential equations describing the probe motion 
is given by r = —VU po t, where r pot represents the acceleration induced by 
Vesta gra vity fie l d and can be expressed as the gradient of the following 
potential ( Kaula . 19661 ) 



U P ot(r, 8,<j>) = 1 — ( — ~) 'Pni^cfyiCnm cos m5+S nm sin mS) , 

(2) 

where the quantities C nm and S nm are the spherical harmonics coefficients 
of the potential of Vesta, n and m are respectively the degree and the order 
of the development, fx and R e are respectively the gravitational constant 
and equatorial radius of Vesta. The quantities (r, 5, eft) describe the plane- 
tocentric spherical coordinates of the space probe. V™ are the associated 
Legendre polynomials. It is worth noting that all terms S n o are zero. The 
potential ([2]) may be expressed using an alternative way, with a cosine terms, 
introducing a phase difference 5 nm as well as a coefficient J nm : 

oo n / ft \ n 

U pot (r,d,(/)) = ( — ) 7^(sin^)(J„ m cos(m5 - m5 nm ) , 

n=2m=0 ^ ' 

(3) 



where 

Cnm — Jnm cos {jTi6 nrJ i)^ S nrn — J nrn s\xi(jYi5 nrr i) ^ TTi5 nrn — arctan 



-s„ 



nm 



Cnrr. 

(4) 

Tricarico and SvkesI (|20ld ) have shown that the dynamical behavior of Dawn 



is dominated by the first 8 degrees of the gravitational potential and that the 
main resonant terms are of degree 3 and 4. Henceforth we model the gravity 
potential of the central body using the coefficients until degree and order 4. 
The coefficients used here can be found in the Table 1 of Tricarico and SvkesI 
(|2010h . 



The numerical integrations of the differential equations are made by a 
home-made numerical software: Numerical Integration of the Motion of 
an Artificial Satellite orbiting a TElluric Planet, for short NIMASTEP. This 
extensive tool allows to derive the osculating motion of an arbitrary object 
orbiting any of the terrestrial planets of our Solar Syste m. This software ha s 



bee n elaborated by N. D elsate and successfully used in IValk et al.l (|2009bl ) 



Compere et al.1 (120111 ) 



and belsate et al.l (l2010h . It has been recently extended for asteroids in 
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3.1 Maps of variation of semi-major axis, inclination and 
radial distance 

In Fig. [2j we plot the results of the numerical integrations performed for a 
set of 316 800 orbits, propagated over a time span of 1 year with a entry-level 
step size of 100 seconds. We consider a set of initial conditions defined by a 
mean anomaly grid of 1° and a semi-major axis grid of 700 m, spanning the 
[385, 1 001] km range. In this work, we explore the dynamics of the spacecraft 
Dawn in a polar orbit (io = 90°) because the initial inclination of Vesta is 
Sun-synchronous, thus corresponds to a quasi-polar orbit. The other fixed 
initial conditions are Q = lo = 0° for the longitude of the ascending node 
and the argument of pericenter respectively, and 9 = 0° for the sideral time. 




100 200 300 100 200 300 o 100 200 300 



Mean anomaly [degree] 

Figure 2: The semi-major axis, the inclination and the radial distance com- 
puted as functions of the initial mean anomaly and the initial semi-major 
axis. The equations of motion include the central body attraction, the har- 
monics until degree and order 4. The anomaly step is 1° and the semi-major 
axis step is 700 m. The initial conditions are iq = 90°, Q = oj = 0° and 
9 = 0°. The integration time is 1 year. The patterns have been obtained 
by plotting the amplitude of variation of the semi-major axis (a amp ) in the 
left panel, the amplitude of variation of the inclination (e am p) in the central 
panel and the amplitude of variation of the radial distance in the right panel 
(ramp)- These amplitude are shown by the colorbars. 

We show the amplitude of the semi-major axis (namely the difference 
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between the maximum and minimum semi-major axis reached during the 
integration), inclination and radial distance of each orbit respectively in 
the left, central and right panels of Fig. [2j These figures allow to locate 
the ground-track resonances. Indeed, when the orbital period of the space 
probe is close to a commensurability with the rotational period of Vesta, 
the repeating ground-track of the probe will periodically encounter the same 
configuration of the gravitational field. This is what we call a ground-track 
resonance. Then a part of the effects of the gravitational perturbations is 
amplified leading to a motion with very long periods and large perturbations 
in the semi-major axis of the spacecraft orbit. This is what we observe in 
the left panel of Fig. [2] where the variation of the semi-major axis is clearly 
largest at ~ 550 km and ~ 720 km, as shown by the color c ode. 

It is obvious that the usual pendulum-like phase spaces (jGedeonl . Il96flh 



(appearing in the case of a ground-track resonance) are distorted like a wave. 
Actually, it is worth noting that, at first, the sampling is carried out with 
respect to the osculating initial condition s. Second, within the fram ework 
of mean motion theory, it is well-known ( Exertier and Metrisl . 19951 ) that, 



due to the short-period oscillations, the mean and the osculating initial con- 
ditions cannot be considered as equal. In other words, for the same fixed 
value of the initial osculating semi-major axis and for various initial mean 
anomalies, we obtain different values for the mean semi-major axis. For 
more explanation in a geostationary satellite case, the reader can refer to 
Valk et al.1 (j2009bh . 



The results obtained here are similar to lTricarico and Svke 2 (|2oioh 



ones 



(Fig. 4 in their article). We easily find 5 main ground-track resonances (noted 
in white in Fig. [2]): 3:2, 4:3, 1:1, 2:3 and 1:2. It is also possible to find other 
resonances (5:4, 4:5, 3:4, 5:3 and 5:2) but they are very thin (noted in gray 
with an arrow in Fig. [2|) . The amplitude of variation of the inclination is near 
1°, except around the separatrix of the 1:1 resonance where the amplitude 
of variation reaches 3.8°. In consequence, in the future, we shall consider 
the inclinations as equal to 90°. In the right panel of Fig. [2j we notice 
that the two most important resonances are the 1:1 and the 2:3 resonances: 
the largest is the 1:1 but the strongest perturbations on the radial distance 
appear in the 2:3 resonance. This last observation will be explained in this 
paper (Section [5]). 
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3.2 Location and aperture of the main ground-track reso- 
nances 

Now let us determine the location and the aperture of the main resonances. 
We perform two initial mean anomaly sections (Fig. [2]) near a stable (vertical 
continued black line) and an unstable (dashed red line) equilibrium point of 
the resonance 1:1. For each orbit of these two sections, we compute, using the 




0.5 0.75 1 1.25 1.5 1.75 2 2.25 2.5 

|Period|x 5.3421288 h 



Figure 3: The mean semi-major axis of a probe around Vesta with respect 
to its mean anomaly period. The periods are obtained by frequency analysis 
(NAFF). The estimation of the periods is made over a period of 1 year. The 
model used is the same as in Fig. [2j The continued black and the dashed 
red lines correspond to the vertical lines of Fig. [2] (left panel). We tick off 
the detectable ground-track resonances. 



Numerical Analysis of Fundamental Frequencies, for short NAFF (jLaskarl . 
20051 ). the fundamental period of the mean anomaly with respect to 



1988, 



the mean initial semi-major axis ao- The mean semi-major axis is calculated 
with a numerical averaging method applied on the osculating semi-major 
axis, on 5 revolutions of Dawn around Vesta. The results are presented in 
Fig. [3] and in Tab. [TJ For more explanation s and other applicatio ns of the 



using of NAFF in space geodesy we refer to iLemaitre et al 
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Table 1: The mean semi-major axis location and the aperture of the main 
resonances (using Fig. [3]). 



Resonance 


3:2 


1:1 


2:3 


1:2 


Mean Semi- major axis [km] 


409.1 


539.6 


717.3 


871.5 


Aperture [km] 


18 


63 


22.5 


7.5 



3.3 Period of the resonant angle at the stable equilibrium of 
the 1:1 resonance 

In the case of a 1:1 ground-track resonance, the resonant ang le is defined by 
a = X — 9 = ("O + U3 + M) — 6 where 6 is the sideral time (jGedeonl . 11969 : 



Valk et all l2009al ). A is the mean longitude, f2 is the ascending node, u) is 



the argument of pericenter and M is the mean anomaly. For illustration, we 
draw in Fig. [J] the motion of a probe in resonance 1:1 with its central body. 






a cos a 



Figure 4: General illustration of the motion of a 1:1 ground-track resonant 
probe near stable equilibria in a rotating reference frame as seen from the 
pole. The variations on the mean semi-major axis have been amplified for 
the illustration. 



We can determine the period of the resonant angle a at the center of 
the banana curves. For that, we use the Namur Algorithm For Forced 
Oscillations, for short NAFFO ( Novelles et al. . 201 ll ). to determine numeri- 
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cally the semi-major axis and the value of the resonant angle corresponding 
to the stable equilibrium. For clarity, the algorithm is reminded here: the 
initial conditions (a, a) are conveniently chosen to be close enough to the 
equilibrium, here for example a = 550 km and a = 90°. Now we try to 
minimize the free libration (the oscillation of a and a in the banana curves) . 
This minimization is performed iteratively by 

(1) performing a numerical integration of the equations of the motion; 

(2) identifying the free component by a frequency analysis (here using 
NAFF); 

(3) removing the free component from the initial conditions to determine 
a new initial condition. Then a new numerical integration is being 
performed (return to step 1). 

We apply this algorithm to different models containing only the gravity field 
of the central body. They respectively contain the J22 term ( J nm are defined 

by Eq. (jJJ)), the J 22 and J2 terms (— C20 = J20 ™= J2), and finally the whole 
gravity field until degree and order 4 (n = 4 = m). When the algorithm 
stops, we make a frequency analysis (NAFF) of the signal (a, a) to determine 
the period of the resonant angle at the equilibrium. 

The results of NAFFO are summarized in Tab. [2j Let us note that we 
calculate an average of the resonant angle and of the semi-major axis with 
respect to the time to obtain a mean value. The 11.21 km difference between 
the J22 and the J22 + J2 model underlines the fact that the J2 term induces 
an important shift on the semi-major axis location. 

Table 2: Numerical study of the location of the 1:1 resonant for different 
models. The first two lines are given by the NAFFO algorithm. The mean 
semi-major axis and the mean resonant angle are obtained by numerical 
average with respect to the time. The third line is obtained using NAFF. 



Model 


J22 


J2 + J22 


All until n = m = 4 


Mean semi-major axis (km) 


549.63 


538.42 


538.89 


Mean resonant angle a (degree) 


90° 


90.814° 


91.144° 


Period (day) 


2.45 


2.47 


2.56 



Our semi- major axis location of the 1:1 s t able e quilibria shows a good 
agreement with Fig. 4 of Tricarico and Sykes ( 2010l ). For the period of the 
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resonant angle, Tricarico and Svkes ( 2010l ) obtained a value close to 2.5 days 
(Fig. 8 in their article) that coincides well with our result. Remark that, for 
this main resonance, our numerical simplified model J% + J22 obtains a good 
approximation of the stable equilibrium location and of the resonant angle 
period. 



4 Ground-track resonances: analytical study 

In this section, we study analytically the main ground-track resonances (1:1, 
1:2, 2:3 and 3:2) found in the previous section. First, we present a general 
Hamiltonian describing the motion of a probe orbiting an irregular central 
body. Then, we identify the main ground-track resonances. For each of 
them we simplify the Hamiltonian in the circular polar case, we explain the 
resonant angle, we determine the location and the aperture of the resonance 
and we calculate the period of the resonant angle at the stable equilibrium. 
The results are compared with the previous numerical results. 
For the 1:1 resonance we also study the circular equatorial case for compari- 
son and we investigate the effect of an error on C20 and C22 on the properties 
of the resonance. 



4.1 Hamiltonian in terms of the orbital elements and condi- 
tion for ground-track resonances 

The potential ([2]) may be expressed using an alternative way, with_the satel- 
lite orbital elements (a, e,i, Q, oj, M) and the sideral time 9, Kaula ( 19661 ) 
represents the gravitational field as 



00 n n +00 



U 



V" \ - \ - \ - \ - 11 ( R e 



+EEEE 



r 1 — ' * — ' < — ' < — ' a \ a 

n=2 m=0 p=0 q=—oo 



F^nmpif) G n pq{&) S nm pq(yi, CO, M, 9) , 

(5) 



where the functions S nmpq depend on the gravity field coefficient C nm and 



Snmpq(Q,U, M,9) 



+ 



+C n m 
Snm 



n—m even 



cos @ nmpq (tl,Lj,M, 9) 



n—m odd 



(6) 



n—m even 



sm& nmpq (n,uj,M, 9) , 



n—m odd 
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where n, m,p, q are integers, 9 is the sideral time and 
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nmpq (n, u, M, 9) = {n- 2p)u + {n-2p + q)M + m(Cl - 9) (7) 



is Kaula's gravitational argument. 

Fnmp(i) and G npq (e) are respectively the inclination and eccentricity 
fu nctions . To ha ve a complete list of the values of these functions, we refer 
to iKaulal (Il966l l and IChaol (2005). The third index q of the eccentricity 



functions determines the lowest power of the eccentricity appearing in the 
potential. 

With the momentum L = ^[JTa conjugated to X = Q + u + M the 
Hamiltonian describing the motion of a probe around an irregular central 
body is given by 



U 



|f+EEE E Rej^Fn mp (i)G npq (e)S nmpq (n,u;,M,9)+9A. 

n=2 m=0 p=0 q=—co 

(8) 



co n n 



n+2 



#A is due to the asteroid rotation, the quantity A is the momentum conju- 
gate to the sidereal time 9. 

The ground-track resonances occur when the rate of Kaula gravitational 
argument is close to zero, that is 



€>nm Pg (A w, M, 9) = (n - 2p)uj + {n-2p + q)M + m{tl 
where p and q are small integer numbers. 







(9) 



4.2 Resonance 1:1 



A similar (more detailed) work has been performed in IValk et al.l ((2009a). 
In this paper the authors studied the geostationary space debris (e = 0, 
i = 0°) without taki ng into account the J2 t e rm. 

It is well know (jOedeoil 1 19691 : IScheeresl . 1 19941 : IValk et all . l2009ah that 
the resonance 1:1 is due to the second degree and order of the gravity field. 
From Eq. (J8|), let us write the Hamiltonian of the second order and degree 
harmonic, denoted by T-Li-.i: 



H 



1:1 



JJ_ 

2L 2 



+ EE E * 

m=0 p=0 q=—oo 



2 M 



e ^ 6 - 2mp 



I ,r.('i) G 2 pq(e) S 2 mpq(tt, u, M, 9) + 9A 

(10) 
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with @2mpq = 2(1 —p)uj + (2 — 2p + q)M + m(0 — 0) . We consider the z-axis 
as aligned with the asteroid main axis of inertia (C 2 i = = S21), q = 
(null eccentricity) and the resonant condition ([9]) becomes: p = 1 
the Hamiltonian writes (S20 = 0) : 



f ; then 



U 



1:1 



2L 2 + e L 6 



^201 00 G 2 w(e) C 2 q cos (9 



2010 J 



(11) 



+ ^22o(i) G 2 oo(e) [C22 cos(6 2 20o) + £22 sin(e 22 oo; 
where 6 nmp9 , F nmp (i) and G npq {e) are presented in Tab. El 



Table 3: Expression of the functions O nmpg , F nmp (i) and G npq {e) with re- 
spect to the values of n = 2, m £ {0, 2}, p = 1 — m/2 and g = 0. 



n 


m 


P 


q 






Gnpq (c) 


2 





1 








| sin 2 i — | 


1 

(l-e 2 )2 


2 


2 








2uj + 2M + 2(fl - 9) 


|(1 + cos i) 2 


(l-§e 2 + ±§e 4 + ...) 



Let us define the so-called resonant angle <r ( Gedeon . 1969 ; Scheeresl . 1994 ; 
Valk et all l2009al ) 



X — 8 = (n + u + M) 



(12) 



In order to keep a canonical set of variables, we use a symplectic transfor- 
mation leading to the new set of canonical variables 



a, 



L' = L, 



A' = A + L. 



(13) 



Then the Hamiltonian (jlip becomes (by selecting only the resonant contri- 
butions or in other words, by averaging over the fast angle 9') 



H 



1:1 



2L' 2 



R 



2 (£ 
L' 6 



I sin 2 i - ~ 



C20 ; 
. (1-e 2 )* 



(14) 



Q r -I I 

+ -(l + cosi) 2 (l--e 2 + — e 4 + ...)(C 22 cos(2 C T) + 5 22 sin(2 C T)) \-0L' . 



A' is now constant and dropped. During the research stages, the orbit 
of Dawn is assumed circular (e = 0). 
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4.2.1 Circular equatorial case: e = 0, i = 

In the circular (e = 0) equatorial (i = 0°) case, the Hamiltonian ([14|) be- 
comes 

Hl:1 = ~2^~ R2e Vs{~l C20+3 ( C22 co < 2 ^ + S ^ sin(2(j)) }-«/. (15) 

Let us note that the eccentricity and the inclination stay constant because 
the Hamiltonian does not explicitly depend on the associated angles u and Cl. 



Scheeres (jl994h studied the stability of bodies in orbits which 



are near- 



synchronous with the asteroid (central body) rotation. If the asteroid is of 
type I, like Vesta, then the stationary orbit contains two stable and two 
unstable equilibria. If the asteroid is of type II, the motion associated with 
near synchronous orbit is always unstable. 

Two stable equilibria!] (cr* l5 L'*^), (cr* 2 , L'* 2 ) as well as two unstable equilibria 
(o"ui,-^ui)> {Puii L'*^) are found to be solutions of 

du = ^ = - (16) 

The values of a* si and a* i are 

a* = iarctanT^ + fc~ k G {0, 1, 2, 3} . (17) 

This formula remains valid for any value of C22 and S22 • For the par- 
ticular value of Vest jg], C m = 3.079667257459264 x 10~ 3 and 5 22 = 
( Tricarico and Svkes . 20ld ). we obtain: 

<t* s1 = k/2, < 2 = 3vr/2 and a* sl = 0, a* s2 = tt . (18) 

The corresponding equilibria in semi-major axis, depending on these values 
and on the value of C20, are the solutions V* of the equation: 

OL' 7 - L'Y ~ \^Rl (C20 ± s ^ + ^'g 2 ' ) = . (19) 

Z \ V G 22 + < b 22 / 

The sign + or — corresponds respectively to the unstable and the stable 
points. The solutions of this equation are given in Tab. HI We point out 



lr The subscrit s for stable and u for unstable. 

2 The corresponding normalized coefficient is equal to C22 = 0.004 771. 



15 



that to neglect the C 2 q term leads to an error spread from 9.5 to 31.9 km 
with a mean value of about 20.7 km on the location of the equilibria i.e. a 
relative variation of ~ 3.76%. For comparison, this difference is about 500 m 
for a GEO i.e. a relative variation of 0.001%. 

Table 4: The circular equatorial case: semi-major axis location of the equi- 
libria (stable a*x = a* 2 "=' a* and unstable a* x = a* 2 a*) of the 1:1 
resonance (second and third columns) with respect to the value of C 2 $ (non- 
normalized coefficient). In fourth and fifth columns the analytical calcula- 
tion of the period (|20p of the resonant angle at the stable equilibrium and 
of the aperture (|22|) of the resonance. 



C20 


a* (km) 


a* (km) 


Period (day) 


Aperture (km) 





544.436 


556.529 


1.216 896 


134.091 


-6.872 554928 x 10" 2 


576.353 


566.066 


1.254 289 


125.726 


Difference 


31.917 


9.537 


134.614/j 


-8.365 



We are now interested in the period of the resonant angle at the stable 
equilibrium (cx*,a*). This is done by linearizing in a neighborhood of the 
equilibrium. Then the Hamiltonian close to the equilibrium is an harmonic 
oscillator that can be expressed in action-angle variables and the frequency 
at the stable equilibrium (the subscript eq. means "evaluated at the equi- 
librium") is given by: 



d 2 n 


d 2 H 


dL' 2 


eq. da 2 



and T = 2tt/&. (20) 

eq. 



The values of the period of the resonant angle at the equilibrium are gathered 
in Tab. HI By a similar approach, we can easily estimate the width of the 
resonant zone; we take the Hamiltonian level curve corresponding to one of 
the unstable equilibria V* and a* noted that we equal to the generic 
Hamiltonian (I15P giving the dynamics a and V along that curve: 

K = -^fp ~ Rljte I - ^20 +3 (C 22 cos(2a) + S 22 sin(2a)) j - 9L' . (21) 

For a = a*, we find the maxima (£ max ) and minima {L 1 ^ ) of this banana- 
curve. We obtain the width of the banana at the stable points: 

AL = -L max — £ m i n , (22) 
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i.e. the aperture of the resonant zone summarized in Tab. 01 
4.2.2 Circular polar case: e = 0, % = tt/2 

The orbit of Dawn is polar and we have seen that the inclination does not 
differ much from 90° for all the radial distances from the center of Vesta 
(Section 13. ip . so let us fix i = it/ 2. 

The methods to find the location of the equilibria, the period at the 
equilibria and the aperture of the resonance are similar to that used for the 
equatorial case (|4.2.ip . 

In the circular polar case, the Hamiltonian (|14p becomes 

= - 2^2 - R l JTs | \ C ™ + \ (°22 cos(2a) + S 22 S m(2a)) j - 0L' . (23) 

The a location of the equlibria is always given by Eq (fT7|) . i.e. for the 
particular value of Vesta 

a* sl = ir/2, <t* s2 = 3tt/2 and a* sl = 0, a* s2 = n. (24) 

The values of the location, period and aperture are gathered in Tab. [5j 

We bring to the fore that, for an inclination of 90°, the period at the 
stable equilibria and the aperture (Tab. [5|) of the resonance are respectively 
equal to the double and the half of the equatorial results (Tab. H]). The 
doubling of the period can be easily understood regarding the equations 
(fl~5j) . (|23|) and (f20j) : in the Hamiltonian (|23j) a \ factor appears, just before 
the C 22 term, (that does not exist in the Hamiltonian f)15|) ) that becomes ^ 
after the square root of the frequency equation ()20|) . which gives a factor 2 
for the period. The aperture can be explained in a similar way. 

We also point out that the location of the resonance is shifted of ~ 7 km 
from the equatorial case for C 2 q = and of ~ 30 km for C 2 q ^ 0. If we 
neglect the C 2 o term, we introduce an error of about 2.6 km on the aperture 
of the resonance and of about 11.8 km on the location of the equilibria. This 
last result coincides well with the 11.2 km numerically found in Tabj2j 

A comparison of the analytical results with the numerical ones shows a 
good agreement (Tab. [H [2] and [5]). The difference between the numerical 
and the analytical results comes from the approximation of the osculating 
elements by the mean keplerian ones. The different apertures of the res- 
onances should also come from contributions that have been neglected or 
averaged in the analytical study. 
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Table 5: The circular polar case: semi-major axis location of the equilibria 
(stable and unstable) of the 1:1 resonance (second and third column) with 
respect to the value of C20 ■ In fourth and fifth columns analytical calculating 
of the period of the resonant angle at the stable equilibrium and of the 
aperture of the resonance respectively. We also recall a part of the results 
of the Tab. [T] and the Tab. for comparison. 



C20 


a* (km) a* (km) 


Period (day) 


Aperture (km) 





552.133 549.113 


2.448 927 


66.699 


-6.872 554 928 x 1(T 2 


540.494 537.159 


2.411514 


69.363 


Difference 


-11.639 -11.957 


-134.686/1 


2.664 


Num. model 


Mean of the s.m.a. 


Period (day) 


Aperture (km) 


<^22 


549.63 


2.45 


66.6 


^22 + C20 


538.42 


2.47 


65.2 


Until degree & order 4 


538.89 


2.56 


63 



4.2.3 Addition of C 40 

Table 6: For the equatorial (noted Eq.) and the polar (noted Pol.) cases: 
semi- major axis location of the stable equilibria of the 1:1 resonance, period 
of the resonant angle at the stable equilibrium and aperture of the resonance 
with respect to the terms used in the Hamiltonian. 



Term 


a* (km) 


Period (day) 


Aperture (km) 


C22 + C20 


Eq. 


566.066 


1.254 289 


125.726 


Pol. 


537.159 


2.411514 


69.363 


C22 + C20 + C40 


Eq. 
Pol. 


567.035 
537.671 


1.253 964 
2.410 739 


125.175 
69.177 



We have seen that the addition of the C20 term has a strong influence 
on the location of the resonance. Then, here, we test the effect of the C40 
term: 

n Jit = H 1:1 - Rj^- Q sin 4 i - f sin 2 i + §) (l + C m . (25) 

We summarize the values of the location, period and aperture in Tab. EJ We 



18 



can see that this term does not significantly affect the results. 



4.2.4 Effect of the error on the coefficient C22 and C20 

It can be shown (B alminol . 1994 ) that 



-*-*>- d§ O'- 2 ^) and c * = m {a '- b2) - (26) 

where a, b and c (with a < b < c) are the semi-axes of the asteroid and R e 
is the mean equator i al ast eroid radius. If we take the values of a, b, c, R e of 
Tricarico and Svkesl lj20ld ): 



a = 289 ± 5km, b = 280 ± 5km, a = 229 ± 5km, i? e = 300km, 

(27) 

then the error on the value of C20 and C22 is 

ACsn ,w A Con AC22 AC22 

— ™ = 18% = and — — = 111% = . (28) 

(^20 C20 ^22 C22 

If we slightly change the reference value i? e or if we add an error (e.g. ±5 km) 
on this reference radius, the results (|28p are similar. The second value of 
([28]) can be surprising but it gives an idea of the boundary on the values 
of C22- Owing to the supposed constant density of the asteroid, the C22 
coefficient is assumed to be positive. Then the minimal (non null) value 
of C22 is arbitrarly chosen equal to 11% of the initial value of C22- So, we 
calculate the extremal values and for each of them we analytically evaluate 
the position (semi-major axis) of the 1:1 resonance (|19j) . the period at the 
equilibrium ([20)) and the aperture of the resonance f[22[) for the polar circular 
case. The results are presented in Tab. [71 The location in semi-major axis 
is between 532 and 541 km. On one hand, for a chosen value of C22, the 
maximum errors on the evaluation of the position of the resonance, of the 
period and of the aperture are respectively equal to 4.6 km, 0.04 d and 1 km, 
due to the error on the C2o- On the other hand, for a fixed value of C20, the 
extremal values of C22 produce errors equal to 3.5 km, 5.6 days and 78 km. 
So it is very important to increase the accuracy on the estimation of the three 
semi-axes (a, b, c). Now if we suppose this improvement, for example an error 
of ±1 km on a, b and c, the error on C20 and C22 becomes respectively 4% 
and 22%. Then the extremal values for the location, period and aperture 
are [536.28 - 538.031] km, [2.181 - 2.733] days and [61.128 - 76.780] km. 
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Table 7: The circular polar case: analytical values of the location (semi- 
major axis noted s.m.a.) of the stable equilibrium of the 1:1 resonance (s. 
s.m.a), of the period at the equilibrium (T) and of the aperture (Apert.) of 
the resonance with respect to the extremal values of the C20 an d C22- The 
middle column and the middle line (T&S) are calculated with the valu es of 
C20 and C22 used in this paper and found in Tricarico and Svkesl (|20ld ) . We 
emphasized the T&S (blue), maximum (red) and minimum (green) values 
obtained for the semi-major axis location, the period and the aperture. 







Min. 


T&S 


Max. 






C20 

C22 ^^--^ 


-0.025 203 


-0.030 735 


-0.036 267 




Min. 


5.248 x 10" 4 


540.870 
7.305 

22.793 


538.655 
7.284 
22.964 


536.397 
7.264 
23.142 


s. s.m.a. (km) 
T(d) 

Apert. (km) 


T&S 


0.004 771 


539.402 
2.418 
68.841 


537.159 
2.412 
69.363 


534.871 
2.405 
69.906 


s. s.m.a. (km) 
T(d) 

Apert. (km) 


Max. 


0.010 067 


537.545 

1.661 
100.221 


535.265 
1.656 
100.991 


532.938 

1.652 
101.194 


s. s.m.a. (km) 
T(d) 

Apert. (km) 



4.3 Resonance 1:2, 2:3 and 3:2 in the polar case 

For the 1:2, 2:3 and 3:2 resonance cases (in the circular polar case), we pro- 
ceed in the same way as in the 1:1 resonance case to find the Hamiltonian, 
the location of the equilibria, the periods at the stable equilibria and the 
aperture of the resonances. 

First, for each of the resonances (1:2, 2:3 and 3:2), we determine the 
values of the integers n,m,p,q to obtain the resonance q\ : q2 to get Kaula 
gravitational argument including the resonant angle a = q2\—q\9. In Tab. [51 
we gather the resonances, the combinations n,m,p,q, the associated reso- 
nant angles, the Kaula gravitational arguments and the functions F nmp {i) 
and G npq {e) respectively evaluated at i = tt/2 and e = 0. 

Second, thanks to this combination of n,m,p,q, we determine the dif- 
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Table 8: Expression of the resonant angles a, the Kaula gravitational argu- 
ments @ n ,m,p,q, the functions F nmp (i = vr/2) and G npq (e = 0) with respect 
to the values of n, m, q = for the resonant case q\ : q 2 - 



Qi ■ Q2 


n 


m 


P 


q 


a 




®nmpq 




G npq {e = 0) 


1:2 


4 


1 


1 





2(0 + uj + M) 


- e 


a-tt 


5/16 


1 


2:3 


3 


2 








3(n + u + M) 


-29 


a-Q 


15/8 


1 


3:2 


4 


3 


1 





2(n + uj + M) 


-39 


a + Q 


105/8 


1 



ferent sets of canonical variables, equivalent to Eq. f)13[) . for each resonance: 

Resonance 1:2 a = 2A - 6, L = 2L' , 9' = 9, A = A' - L'(29) 

Resonance 2:3 a = 3A - 29, L = 3L' , 9' = 9, A = A' - 27(30) 

Resonance 3:2 a = 2X - 39, L = 2L' , 9' = 9, A = A' - 3J$1) 

Then we can write the Hamiltonians assiocated to these resonances: 

u - ^ /x 4 R 2 C 2 o (r{C) s 
Hi* - -^2" 256 jfi (32) 



5 fi 6 R t 



5 4 i cos(a - Q) + C 4 i sin(fj - fi)) + 0(h! - L') 



16384 L' w 

^2:3 - -^72 - 2916L ,6 ( 33 ) 



5p3 



5 

17496 L m 



S 32 cos(ct + C 32 sin(fj - Slf) + 9(A' - 2L') 



H 2 // R\ C 20 , . 

^3:2 - -^2" 256L .6 ( 34 ) 
105 z; 5 /? 4 / \ 

' ' ' S , 43 cos( C 7 + 0) + C43sin( ( j + 0)) + 0(A'-3L')- 



8192 L /10 



In comparison to the Hamiltonian (|23p . the term in front of the cosine is 
—S nm instead of C nm . That is due to the odd value of the n — m (4 — 1 = 3, 
3 — 2 = 1, 4 — 3 = 1) taking place in equation ©. 

With the same procedure as in the 1:1 resonance, we can determine the 
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position of the stable equilibria in the a variable: 



Reso. 1:2 a s = O - arctan ( ) + (2k + l)vr (35) 



5. 



ii 



fc 32 \ 

Reso. 2:3 a s = £1 — arctan I — — I + 2kir (36) 

Reso. 3:2 a s = -Cl - arctan (j^J + ( 2k + l )n ( 37 ) 

and the unstable equilibria: 

Reso. 1:2 a u = VL — arctan ( -5— ) + 2kir (38) 

Reso. 2:3 <r u = Vt - arctan fy^J + 2(A: + l)vr (39) 

Reso. 3:2 <r u = -fi - arctan (^jr^j + 2fc7r • ( 40 ) 

For clarity, we give the mean longitudes corresponding to the stable equi- 
libria (A s ), in degrees: 

a 1 q 

Reso. 1:2 X s = 129.05° + + fc 180°, fc G {0, 1} (41) 

20 4- Q 

Reso. 2:3 X s = 26.45° + + kl20°, fee {0,1, 2} (42) 

3 

on q 

Reso. 3:2 X s = 47.47° + — - — + k 180°, k G {0, 1} (43) 

The values of the locations, periods and apertures of the resonances are 
summarized in Tab. [9l 

We notice that, if we neglect the C20 terms, we introduce a mean error 
of ~ 11 km on the location of the resonances, of ~ 0.5 days on the periods at 
the stable equilibria and of ~ 2 km on the apertures of the resonances. If we 
compare (location and aperture) these analytical results with the numerical 
results (Tab. Q] repeated in the Tab. [9|), we point out that the analytical 
results are adequate. 



5 Increase of the eccentricity in the 2:3 resonance 

In section HI we have analytically found the positions of the main resonances, 
their apertures and their periods at the stable equilibria. As shown in Fig. HI 
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Table 9: The circular polar case: resonances 1:2, 2:3 and 3:2. Values of the 
semi-major axis location (stable and unstable equilibria), the period of the 
resonant angle at the stable equilibrium and the aperture of the resonance 
with respect to the terms used in the Hamiltonian. For comparison, the 
results of the numerical test of the Tab. [T] (noted by Num. in the first 
column) are recalled. 



Reso. 


a* (km) a* (km) 


Period (day) 


Aperture (km) 




C 20 = 


874.090 874.057 


76.341 231 


6.796 


1:2 




866.882 866.847 


75.251862 


9.953 




Num. 


871.5 




7.5 




c 20 = 


721.698 721.365 


9.755 028 


21.953 


2:3 


c 20 /o 


712.898 712.541 


9.608 888 


22.575 




Num. 


717.3 




22.5 




c 20 = o 


420.508 419.911 


4.147596 


20.056 


3:2 


c 20 /o 


404.528 403.755 


3.885 027 


22.401 




Num. 


409.1 
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the ground-track resonances show an oscillation of the semi- major axis and 
of the resonant angle. This is why, in the left panel of Fig. [2j we see an 
increase of the variation of the semi-major axis exactly at the initial semi- 
major axis, corresponding to the ground-track resonances. Up to now we 
have not explained that the strongest perturbations on the radial distance 
come from the 2:3 resonance (right panel of Fig. [2]). In this section, we 
investigate this phenomenon. 

We make a huge number of numerical tests but we present here only the 
significant ones. In Fig. [5j we report the results of a numerical integration 
for a set of 34 920 orbits, propagated over a time span of 1 year. We con- 
sider a set of initial conditions defined by a mean anomaly grid of 1° and a 
semi-major axis grid of 400 m, spanning the [700, 738.8] km range around 
the 2:3 resonance. The other fixed initial conditions are iq = 90° for the 
inclination, Q = to = 0° for the longitude of the ascending node and the 
argument of pericenter, and 6 = 0° for the sideral time. The color code rep- 
resents the maximal value of the eccentricty variation (e max — e m i n ) during 
the integration. 

In the left panel of Fig. [5l we only use the J32 term to modelize the 
central body. As expected by Eq. (j42j) . we see three pendulum-like phase 
spaces. In the central panel, we add the J2 and the J22 terms. We notice that 
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100 200 300 100 200 300 100 200 300 

Mean anomaly [Degree] 



Figure 5: Around the 2:3 ground-track resonance. The amplitude of varia- 
tion of the eccentricity computed as a function of the initial mean anomaly 
and the initial semi-major axis. The anomaly step is 1° and the semi-major 
axis step is 400 m. The initial conditions are in = 90°, = oj = 0° and 
9 = 0°. The integration time is 1 year. The equations of motion include the 
central body attraction, the harmonic J32 in the left panel; the harmonics 
J22 and J32 in the central panel; and all of the harmonics untill degree and 
order 4 (AII44) in the right panel. 
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the eccentricity reaches higher values. In the right panel, the gravity field 
contains all of the terms until degree and order 4. The results are similar 
to these of the central panel. Then, with the help of Fig. [2] and Fig. we 
can conclude that the high value of the amplitude of variation of the radial 
distance that we see in Fig. [2] is due to an increase of the eccentricity in the 
2:3 resonance. We can also deduce from Fig. [5] that the major part of this 
increase of the eccentricity is due to the J22 coefficient. 

Now we analytically justify the fact that the first coefficient of the gravity 
field involving the eccentricity is the J22 term. First, we determine the 
values of the integers n, m,p, q to obtain the 2:3 resonance. Second, we only 
keep the non-null Kaula gravitational argument including the resonant angle 
a = 3X—29. In Tab.fTOlwe present the combinations n, m,p, q, the associated 
Kaula gravitational argument and the functions F nmp (i) (evaluated at i = 
7r/2) and G npq (e). 



Table 10: Expression of the functions Q n mpq 7^ 0, F nmp (i = vr/2) and G npq (e) 
with respect to the values of n, m,p, q for the resonance 2:3. 



11 


m 


P 


q 






G n pq (c) 


2 


2 





1 


3A 


-2B-u- 


ft 


3/4 


7e/2 - 123e 3 /16 + . . . 


3 


2 








3A 


-29- si 




15/8 


1 - 6e 2 + 423e 4 /64 + . . . 


4 


2 





-1 


3A 


-29 + uj- 


SI 


-105/32 


-3e/2 + 75e 3 /16 + . . . 


4 


2 


1 


1 


3A 


-29 -tu- 


SI 


-15/8 


9e/2-3e 3 /16 + ... 



The first line contains the first function G npq directly proportional to 
e, it is related to the J 22 coefficient. The second line corresponds to the 
coefficient J32 that induces the 2:3 resonance. In this term, the eccentricity 
takes place at the second order e 2 . The last two lines are also directly 
proportional to e but, for Vesta, the coefficient C44 is 30 times smaller than 
the coefficient 622- 

If we add the coefficients C22 and S22 to the Hamiltonian (j33|) and if we 
keep only the zeroth and first orders in eccentricity, we obtain: 

^2:3 = -^2-^^(-S 3 2Cos(a-Sl) + C 3 2sm(a-Sl)) (44) 

21u 4 i? 2 / \ 

" ' 1 C22 cos(<7 - oj - Si) + S22 sin((7 - uj - SI) J + 9A . 



8L 6 

With this Hamiltonian we can qualitatively give the secular equation of 

du 



the eccentricity e ~ 1 5^ | . To have an idea of the variation of the eccentricity 
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during the motion we can look at the e/e value: 



&H 



doj 



1 21^R 2 



e 8L 6 



C22 sin(<T — ijj — fi) — S22 cos(<T — ijj — O)^ . (45) 



Then the variation of the eccentricity depends on one of the biggest coef- 
ficients of the gravity field C22- For Vesta, this coefficient C22 is 5 times 
greater than the coefficient C32 responsible for the 3:2 ground-track reso- 
nance. We evaluate (for the values linked to Vesta), near to the 3:2 reso- 
nance, the multiplying coefficient appearing in front of the two major con- 
tributions C32 and C22: 

^5c 32 «6913 and e ^?C 22 « e 25179 . (46) 

The second value is 3.6 times greater than the first. In the case of Dawn 
orbiting Vesta, this highlights the importance of this additional term C22 in 
the 2:3 resonance. 



6 Thrust and probability of capture 

Due to the low thrust propulsion, Dawn will slowly cross the different 
ground-track resonances previously detected. During this transit, the or- 
bit could become chaotic, or Dawn can be trapped in a resonance. 

In this section, we first reproduce the results of Tricarico and Sykesl 



that we explore more in depth. We simulate the slow descent of 



Dawn from HAMO, using a thrust of 20mN, to LAMO. An initial polar 
circular orbit and a continuous thrusting are assumed, with the direction of 
the thrust opposite to that of the relative velocity of Vesta. The simulations 
are displayed in Fig. [6l In the right panel we draw the evolution (continued 
lines) of the (osculating) semi-major axis with respect to the time. The dark 
and light continued curves correspond to an initial mean anomaly respec- 
tively equal to 272° and 250° . On the same panel we plot with dashed lines 
the corresponding evolution of the mean semi-major axis. In the left panel, 
the gray background represents the eyes of the resonance 1:1. To clarify, we 
plot the evolution of the mean semi-major axis in the phase space (a, a). We 
notice that the capture in resonance depends on the initial mean anomaly. 
Actually the probability of capture depends on the phase of the resonant 
angle at the time of the separatrix crossing. The trapping occurs when 
the spacecraft crosses the 1:1 resonance near the longitude where the Vesta 
equatorial semi- axis is the shortest. The reason is that this is for this value 
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Figure 6: Evolution of Dawn from HAMO to LAMO with a thrust of — 20mN 
with two different initial mean anomalies. In the right and left panels, the 
dashed lines (respectively the continued) represent the evolution of the mean 
semi-major axis with respect to the time (respectively to the resonant angle). 
In the right panel, we overlay the evolution of the osculating semi-major axis 
(continued lines). The trapping in the resonance 1:1 depends on the phase of 
the resonant angle at the time of crossing the separatrix. The force model 
contains the gravity field until degree and order 4 and a thrust of 20mN 
opposite to the Vesta rotation. 
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(<7 = 90°) that we analytically find the stable point of the 1:1 resonance 
(Eq. [24")) . In Fig. we plot the orbit (dark line) , which leads to a capture 
in the 1:1 resonance, in a polar coordinate projection. We rediscover the 
banana curve of the Fig. [U 
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the 2:3 resonance, pumping up the eccentricity of the orbit. The right panel 
of our Fig. [2] also presents these structures and we have explained the origin 
of this phenomenon in the previous section. So it is possible that the increase 
of the eccentricity affects the probability of capture. Let us calculate the 
probability of capture with respect to the thrust and the initial semi-major 
axis. The results are presented in the Fig. [8l To obtain these figures, we 
take a force model that contains the gravitiy field until degree and order 4 
and a thrust. The thrust is spread from 20mN to 36mN (to 50 for the right 
panel) with a step of 0.5mN. For each value of the thrust, we numerically 
integrate 50 times 1 000 orbits with a random initial mean anomaly. In the 
left panel, the initial semi-major axis is equal to 1 000 km corresponding to 
the HAMO altitude. Then these orbits cross the 2:3 resonance. In the right 
panel the initial semi-major axis is equal to 600 km and does not cross any 
resonance except the 1:1. 
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Figure 8: Probability of trapping of Dawn in the 1:1 resonance with respect 
to the thrust and the initial semi-major axis. For each value of the thrust 
we make numerical integrations with 1 000 different mean anomalies chosen 
randomly. We repeat this operation 50 times to obtain the means and the 
standard deviations. 



In the left panel, we notice that the probability of capture is not linked 
to the value of the thrust. The mean pro bability of capture i s equa l to 
8.26% ± 3.288. This confirms the results of lTricarico and Svked (|2010l ). In 



the right panel, the probability of capture fluctuates between 4% and 18% 
with a mean equal to 10.24% and a standard deviation of 3.688. We observe 
that the probability of trapping is directly linked to the value of the thrust. 
From a thrust equal to 48mN, the probability decreases. So, for any initial 
semi-major axis and mean anomaly, we advice to take a thrust between 20 
and 31mN or greater than 48mN. In this way, we obtain an higher proba- 
bility to cross, without permanent trapping, the 1:1 resonance. 

Now, if the probe is trapped in the 1:1 resonance, it is interesting to 
know how we can get the probe out of the resonance. First, in Fig. [6l 
for the captured orbit (dark line) in the eye of the resonance, we carry on 
the thrust of 20mN. Thus, to escape from this resonance, we increase the 
val ue of the thrust. We dec ide to increase the thrust from 20 to 35mN 
(as Tricarico and Sykesl 2010l ). We perform two tests where we modify the 



thrust at two different times: after 12 and 13 days. We note that the orbit 
with the change of thrust at 13 days stays in the resonance 1:1 and is closer 
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Figure 9: From a capture in the 1:1 resonance (black line), thrust of 35mN 
at two different times to try to escape from the trapping. The force model 
contains a thrust and the gravitiy field until degree and order 4. In right 
and left panels the dashed lines (respectively the continued) represent the 
evolution of the mean semi-major axis respectively with respect to the time 
and to the resonant angle. In right panel we overlay the evolution of the 
osculating semi-major axis. The black line correspond to the dark line of 
the Fig. [H The escape from the resonance 1:1 depends on the phase of the 
resonant angle at the time of the operating phase (changing of the thrust). 
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to the center of libration. The other case escapes from the resonance. The 
difference in the initial time of the 35mN thrust induces a different value of 
the phase of the resonant angle at the time of the change of the thrust. These 
maps (the left panels of Figs. [9] and [6]) allow us to easy understand that the 
probability of trapping or escape from the 1:1 resonance is dependent on the 
resonant angle phase. 



7 Conclusions 

The orbit dynamics of a polar space probe orbiting Vesta has been investi- 
gated. The proposed model includes the effects of the spherical harmonic ap- 
proximation of the gravity field of the asteroid. We study the main ground- 
track resonances: 3:2, 4:3, 1:1, 2:3 and 1:2. We have developed a general 
method using an averaged Hamiltonian formulation. For each main ground- 
track resonance, we have analytically localized the resonance, determined 
the aperture and calculated the period at the stable equilibria. We have 
also discussed the effect of the error on the C20 and C22 coefficients on the 
properties of the 1:1 resonance. 

Our method allows to have an analytical global approach to search all 
main ground-track resonances and to compute the location and the aperture 
of the resonances. The analytical results have been checked and validated 
numerically by performing numerical integrations of the complete systems. 

Our theory is a b le to reproduce, explain and complete the results of 
Tricarico and Syke 1 (l201(]h . The theory is general enough to be applied to 



a wide range of probes around any dwarf planet or most massive asteroid. 
In opposition to a numerical study, if the values of the spherical harmonic 
coefficients change (for example after new observations) our method stays 
valid. 

In a numerical way, we have found and justify that the increase of the 
eccentricity in the 2:3 resonance is due to the C22 coefficient. We have also 
numerically studied the probability of trapping in the 1:1 resonance that 
does not depend on the thrust if the initial semi-major axis of Dawn is 
upper than the 2:3 resonance location. We have also showed that the phase 
of the resonant angle is the key during the operating phase. 
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